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Abstract 

We study the sensitivity of an asymmetric i3-factory to indirect CPT violat- 
ing effects. We find that both dilepton signals and nonleptonic asymmetries 
can be used to bound the space of parameters describing possible CPT vio- 
lation in the mixing matrix of the neutral B system. 
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I. INTRODUCTION 



Invariance under CPT transforms is an exact property of a relativistic, local and renor- 
malizable field theory, satisfying causality and invariance under proper Lorentz transforma- 
tions, incorporating usual spin-statistic relations and provided with asymptotic states [[!]. 
Therefore, tests of CPT invariance probe fundamental aspects of the current field-theoretical 
description of microscopic phenomena, and provide stringent constraints on quantum grav- 
ity, supergravity and superstring-inspired scenarios requiring a violation of CPT 0. 

CPT symmetry, which implies the equality of mass, lifetime, branching ratios of particles 
and antiparticles, was tested at various levels of accuracy since the early days of the high- 
energy physics era |J. The main activity was focused on the theoretical and experimental 
analysis of the neutral kaon system, where it was observed that the phases 0oo an d 
4>sw are equal within an accuracy of 1° 0]. At ^-factories, where 8.9 x 10 9 — K$ pairs 
per year will be produced, tests of even higher accuracy are expected to be performed in the 
near future 0. 

Tests of CPT invariance have been proposed for the neutral D meson system ||. More- 
over, the symmetric and asymmetric 5-factories currently under construction |7|-|| represent 
experimental facilities where CPT symmetry can be investigated in the neutral B system. 
In this case, information is used from T(4S) decays to coherent C-odd B°-B° states (C-even 
states, e.g. from T(4S) — > £> _B°7, are expected at the level of 1/10 9 events [0) so that 
clear CPT tests, in addition to CP measurements, can be carried out. 

The proposal of using the neutral B system for testing CPT was already put forward in 



refs. [|14|12|j . In particular, in ref . |L2| the features of symmetric and asymmetric I?-factories 
were considered in order to put bounds to possible violations of CPT, mainly using the 
channel where one neutral B meson decays via a semileptonic transition, and the other one 
decays to J/ipKs; Monte Carlo simulations were used to study the sensitivity to this process 
in various experimental environments. In this paper we want to reconsider the problem. In 
particular, we investigate the feasibility of bounding CPT invariance in the neutral B system 
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at an asymmetric B-factory using not only nonleptonic decay channels, but also the decay 
modes where both the neutral B mesons decay via a semileptonic transition, since a high 
statistics is expected in this case, with limited reconstruction and background rejection 
difficulties. We introduce a CPT violating term in the mixing matrix of neutral B mesons, 
according to the notation in ref. JTTJ, and then we consider the case where lepton pairs are 
identified in semileptonic B° and B° decays. 

We also consider the decay channel B — > J / ipKg and report on the sensitivity to different 
decay modes, such as B — > D + D~ and B — > D* + D*~ . 

Our conclusion is that it is possible to sensibly constrain the parameter space of CPT 
violating effects in the neutral B mixing matrix. The future 5-factories, therefore, represent 
powerful facilities for probing a fundamental aspect of the field theoretical description of the 
elementary interactions. 

II. PARAMETERIZING CPT VIOLATION IN THE B° - B° MIXING MATRIX 



In the Wigner-Weisskopf approach, the mixing B°-B° is governed by the Hamiltonian 

H = M-i-=[ * 2 (2.1) 

2 \M* 2 — M 22 -i v -f) 

f B °\ - - 

acting on the vector . The condition H n ^ H 22 , or < B°\H\B° > ^ < B°\H\B° > , 

\B° 



implies an indirect violation of CPT invariance in the neutral B system. Following |IT] and 



|13|| , such a violation can be easily parameterized expressing the matrix (|2.1|) as follows: 

H = —iD + o-i^! + a 2 E 2 + a 3 E 3 (2.2) 

with (jj the Pauli matrices and D, E%, E 2 , E 3 four complex numbers. The correspondence 
between the parameters in eqs.(2.1- 2~2]) is simply given by: 



E x = ReM 12 - % -ReY 12 
E 2 = -ImM 12 + % -ImT l2 

E 3 = i(Mn - M 22 ) - % -{Y 11 - T 22 ) . (2.3) 



CP and CPT violating effects can be parameterized []TT| , p!3| by introducing the new complex 
variables E, 9, 4>, defined as: 

E = (E\ + E% + El) 1 ^ 2 , £1 = sin 6> cos , E 2 = E sin 9sm(p , E 3 = E cos 9 . (2.4) 

As a matter of fact, CPT symmetry implies M n = M 2 2 and r n = r 22 , i.e. cos# = 0. 
Moreover, using the phase convention CP\B° >= \B° >, CP invariance requires the equality 
< B°\H\B° >=< B°\H\B° >, i.e.: 

M 12 -* T -f = M* 2 - i^f . (2.5) 

This means Im Ti 2 = Im M\ 2 = 0, i.e. = and cos 6* = 0. 

The coefficients relating the mass eigenstates of the Hamiltonian (|2.1| ) to the flavour 

(-) 

eigenstates | B°>: 

lft >= w? + wy (p2|s ° > ~" 2}B ° >) (27) 

are simply related to the parameters 9 and <ft by the equations 

£ = e *tan? . ^ = e*cot-. (2.8) 
Pi 2 p 2 2 1 ; 

Then, CPT invariance implies # = | and 

* = » = e** ; (2.9) 

Pi £>2 

a violation of CPT is governed by the parameter 

s = cot9 = -( q ^- q Ae-^ . (2.10) 
2U Pi/ 1 ' 

Although, in principle, 4> an d 9 are complex numbers, we assume in the following 

Im<j) = , Im 9 = , (2.11) 



and therefore a real variable s, to reduce the phenomenological analysis to the effects of the 
variation of a single parameter. We shall comment below on the possible effects related to 
the imaginary part of s. 

If at t = only a B° or B° state is present, the time evolution of the mass eigenstates 
l-Bi > and \B 2 >: 



\B n (t) >= e-^-^lB^O) > (n = 1, 2) 



(2.12) 



implies that 



B°(t) >= e" (fm+ ^ ^(t)!^ > +g+{t)\B° > 



(2.13) 



B { \t) >= e - (im+ ^[p_(t)|5 u > +g-{t)\B {) >] 



(2.14) 



where 1 1 



g±(t) = cos 2 °- e±(*Am-4F)| + sin 2 ? e T(^m-M )f 



(2.15) 



g±(t) = sin -cos - [e (iAm a ) a 



rai(ri) and m 2 (r 2 ) are the mass (width) of \Bi > and \B 2 >, respectively, and m 



(2.16) 



mi+m2 
2 ' 



Am = m 2 - mi, T = AT = I\ - T 2 . 

Since 10" 2 [|TJ, in the following we neglect the actual value of AI\ 

Eqs.( 2.13 - |2.16 ) show that the time evolution of neutral B states is governed also by the 
CPT violating parameter s, and then it is possible to define observables sensitive to s. A 
bound on s can be inferred from the mass difference B°- B°, using the relation 



m#o — m^o 



:(mi - m 2 ) 



(2.17) 



with mi — m 2 = 3 x 10 GeV H the measured mass difference between B\ and B 2 . A 



value s ~ O(10 *) implies a relative mass difference 



\m B o-m B o 



kaon system the bound |t " x ° m,?o1 < 9 x 10~ 19 

average 



average 

corresponds to s < 2 x 10 -4 . 



C(10" 15 ); in the neutral 



III. DILEPTONS 



Let us now consider the signals of a possible CPT violation in various final states, starting 
from the processes where both the neutral B mesons decay via a semileptonic transition. 
The B°B° wave function from T(4S f ) decays reads 



\ip(t) >= -= \\B°(k,t) > \B°(-k,t) > +C\B°(-k,t) > \B°(k,t) > 
V2 L 



(3.1) 



— * 

where ±k are the meson momenta in the T(AS) rest frame, and C = — 1 is the charge 
conjugation of the B°B° pair. 

Let us determine the probability of opposite sign dilepton production (l + X~~), (l~X + ), 
X* 1 being a generic hadronic state. If t\ and t 2 are the decay times of the B mesons 
with momentum k and —k, respectively, and if P ± ^(ti, ^2) is the probability for a positive 
(negative) lepton at the time tifo), one has: 



p- + (t 1 ,t 2 ) oc I < rX + ,t 1 ;l + X-,t 2 \^ > | 2 

< rX + ,t 1 \H w \B°(t 1 ) >< l + X-,t 2 \H w \B°(t 2 ) > 

2 

- < rX+^^HwlB ^) >< l + X',t 2 \H w \B°(t 2 ) > 



(3.2) 



1 

^2 



Since, within the Standard Model, only the transitions B° — > / X + , B° — > l + X occur at 
the leading order in the Fermi constant, the following relations can be derived: 



< l + X-,t\H w \B°(t) > = e ~ (im+ ^5 g_(t) < l + X-\H w \B° > , 

< rX + ,t\H w \B°(t) > = e - (im+ ^ g + (t) < rX + \H w \B° > , 

< rx + ,t\H w \B°(t) > = e -( <ro +5)l g _(t) < rx + \H w \B° > , 

< l + X-,t\H w \B°(t) > = e -( <ro +5)l g+ (t) < l + X~\H w \B° > . 
Neglecting violation of CPT in the decay amplitude, we can put: 

< l + X-\H w \B° >=< rX + \H w \B° >= A, . 



(3.3) 
(3.4) 
(3.5) 
(3.6) 



(3.7) 



Then, keeping only terms of 0(s 2 ), the probability of production of opposite sign dileptons 
is given by: 
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P- + (h,t 2 ) cc \e- T ^\\A l \ 2 f\Mti)9-{t2) - 9-{ti)9 + {t 2 )\ 2 

I e -r(ti+t 2 ) 
= -(|^| 2 ) 2 — — [1 + 2s 2 + cos Am{t x - t 2 )} . 



(3.8) 



1 + s 2 

The production probability of equal sign dileptons can be derived in an analogous way; for 
C = — 1 one has 

-r(t 1 +t 2 ) 



P^t^ocd^! 2 ) 2 



4(1 + s 2 



■[1 — cos Am(ti — t 2 )] 



(3.9) 



For s = 0, eqs. (|3.8| - |3.9|) reduce to well known expressions for the production probability of 
lepton pairs in B°-B° decays. Time integration of (|3.8|) gives the integrated probability 



oc 



1 2\2 



1 



where Xd 



Am 



P 2(1 + s 2 ) 



The quantity ^ A ^[ ^ is proportional to the product of the semileptonic 



1 + 2^ 



1 



1 + x 2 



(3.10) 



B°, B° branching ratio; the fractions n ±± of equal and opposite sign dileptons produced for 
each decay can be easily derived: 



n 



2(1 + s 2 ) 



1 + 2s 2 + 



and 



n 



n 



l + x 2 . 



(3.11) 



4(1 + s 2 ) I + x 2 d ' 
The dependence on the parameter s signals CPT violation. 
A measurement of the ratio 



(3.12) 



R 



n 



n 



IV 



(3.13) 



l + (l + 2s 2 )(l + x 2 ) 

provides a bound on s. Moreover, a measurement of time-dependent production, possible at 
an asymmetric .B-factory, provides us with the fraction of opposite and equal sign dileptons: 



n~ + (At) 



-rlAil 



2(1 + s 2 ) 



[1 + 2s 2 + cos AmAt] 



n ++ (At) = n~~(At) 



-r\At\ 



4(1 + s 2 
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■[1 — cos AmAt] 



(3.14) 



(3.15) 



Then, the ratio R(At), analogous to ( |3.13|) , can be written as: 

= n ++ (At) + n-(At) = 1- cos Am At 
1 ' ' n-+(At) " 1 + 2s 2 + cosAmAt ' l ' ' 

It is worth observing that other asymmetries, obtained reconstructing the temporal se- 
quences of production of equal and opposite-sign dileptons: 

n ++ — n — n~ + — n H 

A=^— , A — , 3.17 

n ++ + n n + + n + 

vanish (both for positive and negative charge conjugation of the initial B° B° states). How- 
ever, this is only the case for a real s. In general, one finds [ljj] : 

= n-+(At)-n+-(At) = 2(1 m s) sinAm|At| 
1 J n-+(At) +n+-(At) 1 + cosAmAt ' l ' ' 



and therefore from the asymmetry ( |3.18|) a bound on Im s can be derived. 



It is now interesting to determine the upper bound on s that is possible to obtain in 
the BaBar experiment at the SLAC B-factory PEP II making use of dilepton events. At 
PEP II N(B°B°) ~ 1.8 x 10 7 B°B° pairs are expected per year, for a machine running at 
the design luminosity Using the value of the semileptonic branching fraction B(B° — > 
l + X~) ~ 0.105, and the estimated lepton tagging efficiency e tag = 0.65 J7|, the number of 
expected dilepton events per year is: 

N ++ = N = N(B°B°)[B(B -> l + X-)] 2 el g n ++ ~ 2.8 x 10 4 , (3.19) 

jV~+ = N(B°B°)[B(B° -> Z+X-)] 2 e 2 as n- + ~ 2.8 x 10 5 . (3.20) 

The error on the variable R can be derived assuming a binomial distribution for A^ ±± and 
A^~ + . If p and 1 — p are the probabilities of two leptons produced with opposite and equal 
charges: 

P = n~ + = 1 (l + 2s 2 + -^ 1 ) , (3.21) 
2(1 + s 2 ) V 1 + 



l- p = n ++ + n — = X X \ / 322) 

2 1 + s 2 l + i 



(for s 2 « 1), one has that the probability distribution 

P(N- + )=( N ) p N ~ + (1 - p) N - N ~ + , (3.23) 
with N = N ++ + A^ + N~ + , can be expressed in terms of R in the form 



N 



N 



1 + RJ \ jv_ 

l+R 



p T+R (1 -p)"-T+5 ; (3.24) 



yielding 



N , 

= Np(l-p) (3.25) 



l + R 



and finally the error a(R): 



9 N~+ N++ 

*(R) = (l + Ry ] J . (3.26) 

In Fig.l the ratio R in eq.(3.13) is plotted as a function of Xd and s in the ranges 
0.63 < x d < 0.83 and < s < 0.3. In Fig. 2 we plot the time dependent ratio R(At) 
in eq.(3.14). The corresponding error bars have been calculated using (3.26). For time- 
integrated measurements, the error a(R) is given by a(R) = 9 x 10~ 4 for Xd = 0.73, while, 
considering the whole range of x d) a(R) runs from 8 x 10~ 4 to 1 x 10~ 3 . 

These errors give an insight on the possible accuracy that can be obtained at a £>-factory 
such as PEP II, due both to the machine luminosity and to the efficient lepton identification. 
The small value of cr(R) shows that the region of possible values of x d and s can be tightly 
constrained. The measurements of time integrated and time-dependent dilepton production 
fractions allow to establish an upper bound to the parameter s. As a matter of fact, from 
Fig. 2 we get that, for AmAt = |, the bounds s < (5, 7, 8) x 10~ 2 are obtained, respectively, 
within one, two and three standard deviations on R(At); for AmAt = | and the bounds 
are s < (4, 6, 7) x 10~ 2 , and s < (4, 5, 6) x 10~ 2 , respectively. 



IV. LEPTON-HADRON CHANNELS 
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A. Decay mode B — ► J/ipKs 



Now we wish to investigate the possibility of bounding CPT violating effects by means 
of nonleptonic decays of neutral B mesons ]E| , using our parameterization of the CPT 
breaking term. We consider events where one B decays via a semileptonic transition to 
Z ± X =F , while the other B decays into a nonleptonic state. On the experimental side, the 
semileptonic decay can be used to tag the flavour of the other B meson decaying into the 
hadronic final state \f >. 

Integrated and time-dependent CP asymmetries of the type 

N(l-X+,f)-N(l+X-,f) 

N(l-X+,f) + N(l+X-J) { ■ ] 

A( At) = W-X + J;*t)-N(l+X-J;At) 
1 1 N(l~X+J;At)+N(l+X-J;At) 1 " 1 

can be used to test CPT invariance. Let us consider, for example, the hadronic decay 

H ... 

B°— > J/ipKs- A tree and a penguin diagram contribute to such a process, so that we can 
express the transition amplitudes in the form: 

< J/^K S \H W \B° >= A T e i&T e i(t>T + A P e iSp e^ p (4.3) 

< J/i)K s \H w \B° >= A T e iST e- i4>T + A P e i5p e-^ p , (4.4) 

where A?,p and <Pt,p are the weak amplitudes and phases respectively, while St,p are the 
strong phases. At the leading order in s the CP asymmetries ( f4.1|) and ( |4.2j ) read respectively 

x 2 

A= cos 2/3, (4.5) 

1 + x* 

A(At) = - sin 2/3 sin Am At + s cos 2/3(1 - cos Am At) , (4.6) 

in terms of the /3 angle of the unitarity triangle. 

The sensitivity of PEP II to a CPT violation effect in this channel can be obtained 
noticing that the error on A is given by 
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a(A) = d, 



+ (47 ) 



eff 



In eg. ( fl7l) the effective number of events N e ff is given by: N e ff = N(f)K(f, f)E tag , with 

N(f) the number of events with a reconstructed hadronic state, K(f,f) = n(f,l + X~) + 

(-) 

n(f, I X + ), n{ f , I X T ) being the fractions of lepton-hadron events produced in the decay 

of the B°B° pair, E tag the total efficiency of tagging, d the dilution factor for the production 

(-) 

of |/ > states from events different from B°— > f decays. For the final state / = J/tpKs 
the values E tag = 0.34, N(f) = 886 and d — 1 are expected 0. In Fig. 3 and Fig.4 we plot 
integrated asymmetries for (3 = 18°, (3 = 22° and j3 = 26°, and time-dependent asymmetries 
for {3 = 22°, with the corresponding error bars. We notice that the maximum value of the 
time-dependent asymmetry grows with s, and therefore higher order terms in the s-expansion 
are needed to fulfil the bound \A\ < 1. 

For integrated measurements at PEP II, the expected error on A is: a = 5.8 x 10 -2 , 
quite independent of (3. As a result, it is possible to bound s < 0.21 for j3 = 18°, s < 0.24 
for f3 = 22° and s < 0.28 for (3 = 26° within one standard deviation. 

From time-dependent asymmetries, in the PEP II experimental environment the bounds 
s < 0.13 (0.26) (within one and two standard deviations) can be established in correspon- 
dence to AmAt = f , and s < 0.17 for AmAi = 2* 



B. Other decay modes 

Other hadronic decay modes, such as B — > D + D~ and B — > D* + D*~ , can be used to 
bound CPT invariance. To calculate the asymmetries for such final states, new information 
has to be considered concerning the transition amplitude and the final state interaction. 

The transition amplitudes B — > £)(*)£)(*) are governed by the effective hamiltonian [15 



H eff (AB = -1) = ^=[V ub V: q { Cl Ol + c 2 O u 2 ) + V cb V: q { Cl Ol + c 2 2 c ) - VfeV^cA] , (4-8) 
where Vjb and V* q (j — u,c,t; q = d, s) are elements of the Cabibbo-Kobayashi-Maskawa 
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matrix, and q (i = 1, ...,6) the Wilson coefficients at the scale fj, ~ m&; 0"' c and 0^ c are 
current- current operators, Oj (z = 3, ..,6) are penguin operators. 

In the vacuum saturation approximation, the decay amplitudes can be written as 

2a 6 m|, 



< D + D-\H W \B° > = —= 



VcdV^ - V td V* a 4 + 



(m b - m c )(m c + m d ) 



x < D+l&f^l - l5 )d\0 >< D-\bY(l - l 5 )b > \B U > , (4.9) 

< D* + D*-\H W \B° > = ^=(V cd V; bai - V td V*a,) 

x < D* + |c T ^(l - l5 )d\0 >< D*~\hHX ~ 7 5 )c|£? > , (4.10) 

with parameters a lr ..,ae related to the Wilson coefficients ci,... ,c 6 by the equations: 

fl2i~i = c 2 i-i + ^ , a 2 j = c 2i + (4-11) 

(i = 1, 2, 3 and A^ c is the number of colours). 

In the NDR renormalization scheme at the next-to-leading order, using N c = 3 and the 



scale /i = mb, we get Jl5| 



ai = 1.017 , a 2 = 0.175 , a 3 = 0.0013 , a 4 = -0.030, a 5 = -0.0037, a 6 = -0.038. (4.12) 

Integrated and time-dependent asymmetries depend on the parameters of the Cabibbo- 
Kobayashi-Maskawa matrix and, unlike the B° — ► J/ipKs mode, also on strong final state 
phases. Various estimates, however, suggest that for B transitions to charmed states such 
strong phases are small |16| , and therefore we neglect them in the analysis. The asymmetry 



( |4.2| ) can be computed in a straightforward way; we omit here the lengthy expression, and 
only plot in Fig.5 and Fig. 6 the time dependent asymmetries for B — > D + D~ and B — > 
D*+D*- respectively, for (3 = 22° and the parameters of the CKM matrix p = 0.05, rj = 0.39, 
7 = 83°; the corresponding error bars are also depicted in the figures. 

The analysis shows that a sensible bound on s can be obtained from these decay channels 
only if tagging and reconstruction efficiencies at PEP II improve with respect to the estimates 
quoted in @ . For instance, an improvement by a factor of two will allow to put the bounds 
s < 0.18 (from B — > DD) and s < 0.22 (from B — > D*D*) in correspondence to Am At = |. 
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As for the hadronic channels B° — > 7r + 7r~ and B° — > p ± 7T T , the CP asymmetries depend 
on the strong phases due to the hadron rescattering in the final state. In any case, we found 
that the statistics expected for such decay modes does not allow to sensitively bound the 
CPT violating parameter s. 



V. TESTING DIRECT CPT VIOLATION EFFECTS 

So far we analyzed indirect CPT violation effects. It is also possible to consider observ- 
ables which are sensitive to direct violation effects. One of them is the asymmetry between 
the partial widths of charged B mesons, related to a CPT violating term in the transition 
amplitudes < f\H w \B + >^< f\H w \B~ >: 

r(g+ -> /) - r(B- -> /) 

r(B+->/) + r(B--/) ' 1 ' } 

We do not parameterize this kind of violation, but just compute the expected error on the 
asymmetry ( |5.1| ). For the decay mode B + — > D°/ + z//, N = 2.9 x 10 5 identified events 



are expected per year 0, so that the estimated error on (O) is a = 7.9 x 10 -3 ; for 
B + — ► D°*l + ui, the number of expected events per year is N = 3.9 x 10 6 , and the estimated 
error is o = 5.1 x 10 -3 . 



VI. CONCLUSIONS 

We investigated the possibility of testing CPT symmetry in the neutral B meson system, 
considering the experimental environment provided by the S-factory PEP II. We found that 
dileptons and lepton-hadron events, with the hadronic decay B° — ► J /ipK Sl can be used to 
constrain the CPT violating parameter s. 

In particular, from equal and opposite sign dilepton production, one can derive the bound 



As for the hadronic decay B° — > J/ipKg, the sensitivity of PEP II is up to s ~ O(10 _1 ) 



s < (9(10 2 ), which corresponds to a relative mass difference ^~ m m\ < O ( 10 



Our analysis, though being quite different, confirms the estimate in |12| about the order of 
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CPT violating effects that can be constrained at PEP II. 

Other channels, such as B° — > D + D~ and B° — > D* + D*~, might allow a sensible 
bounding of s if tagging and reconstruction efficiencies improve at least by a factor of two 
with respect to the current expectations. Charmless decays of the type B° — > it + it~, B° — > 
p ± 7r T , turn out to be unsuitable for testing CPT. 
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APPENDIX A: CASE C = +1 



For the sake of completeness, we report the relevant formulae corresponding to a C even 
B°-B° state: 



I e -r(ti+t 2 ) 
P; + oc i (|A| 2 ) 2 4(1 + s2)2 {l + cos Amfa + f 2 ) 

+ s 2 [l + 2 cos Amti + 2 cos Amt 2 — cos Am(ti — t 2 )]} ■ 



(Al) 



P^ihM) oc (|A| 2 ) 2 ^^{1 - cos Amfa + i 2 ) 

+ s 2 [3 — 2 cos Amti — 2 cos Amt 2 + cos Am(ti + t 2 )] } 



(A2) 



2(1 + s 2 y 



xj + xj + 2 , 2 4 + x 2 



(1 + ^) 2 



1 + x 2 , 



(A3) 



+ 



4(1 + s 2 ) 



™4 I q™2 q™2 



(i + x d y • ~ i + x d 



(A4) 



n~ + (At) 



e- r l Af l f 2 + x 
2(l + s 2 ) 2 ll + a; 2 

+ cos(Am|At| - a) 



d +s 2 



1 — cos AmAt + 



1 + a; 2 /4 



(sin(Am|At| - a) 



(A5) 



n+ + (At) = 



e -r|At| 



4(l + s 2 ) 2 ll+^ 
+ cos(Am|At| - a)) 



3 + cos AmAt + 



V 1 + *2/4 



(sin(Am|At| - a) 



(A6) 



with ck = arctan ^. 
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FIGURE CAPTIONS 



Fig. 1 

The ratio R in ( |3.13| ) versus Xd and s. The continuous line corresponds to s = 0, the 
dashed-dotted line to s = 0.1, the dashed line to s — 0.2, the dotted line to s = 0.3. 

Fig. 2 

Time dependent results for R (eq.( |3.16| )). Notations as in fig.l. The expected statistical 



error is smaller than the size of the dots for the different values of AmAt. 



Fig. 3 



Integrated asymmetry (4.5) for the channel B° — > J/ipKs- The three lines correspond to 
the values of the (3 angle: (3 = 18° (dashed-dotted line), (3 = 22° (solid line) and (3 = 26° 
(dashed line). Statistical errors show the minimum value of s accessible from this mode. 

Fig. 4 

Time dependent asymmetry (4.6) for the channel B° — > J/i[)Kg. The lines correspond to 



the values s = (solid line), s — 0.1 (dashed-dotted line), s = 0.2 (dashed line) and s = 0.3 
(dotted line). The expected statistical error for different values of AmAt is also shown. 

Fig. 5 

Time dependent asymmetry ( f4.2|) for the channel B° — ► D + D~ . The lines correspond to 



the values s = (solid line), s — 0.1 (dashed-dotted line) and s = 0.2 (dashed line). 
Fig. 6 

Time dependent asymmetry (|4.2|) for the channel B° — ► D* + D*~ . The lines correspond to 
the values s = (solid line), s = 0.1 (dashed-dotted line), s = 0.2 (dashed line) and s = 0.3 
(dotted line). 
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